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^ ■ Abstract 

I ■ A modified Kardar-Parisi-Zhang (KPZ) equation is introduced, and solved 

exactly in the infinite-range limit. In the low-noise limit the system exhibits a 

^ . weak-to-strong coupling transition, rounded for non-zero noise, as a function 

. of the KPZ non-linearity. The strong-coupling regime is characterised by a 

. double-peaked height distribution in the stationary state. The nonstationary 

. dynamics is quite different from that of the stationary state. 
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The field of kinetic rougliening of surfaces has seen an explosion of interest in recent years 
(see, e.g, [^0). Numerous theoretical models have been proposed, of which the simplest 
non-trivial example is the Kardar-Parisi- Zhang (KPZ) equation 



where h{x, t) specifies the height of the surface at point x at time t, and r]{x, t) is a zero mean 
Gaussian noise with correlator {ri{x,t)ri{x,t)) = 2Dq6'^{x — x)6(t — t). This same equation 
is believed to describe the statistics of directed polymers in a random medium where 
h has the meaning of a free energy. 

For spatial dimensionality d > 2, the KPZ equation is known to exhibit, for fixed noise 
strength Dq, a phase transition at A = Ac from a 'weak-coupling' or 'smooth' phase for A < Ac 
to a 'strong-coupling' or 'rough' phase for A > Ac- For d < 2, only the strong-coupling phase 
exists for any A > 0. 

Despite intense effort in recent years, the properties of the strong-coupling phase are 
rather poorly understood. In the present Letter, we attempt to rectify this by presenting an 
exact solution of the model in the infinite range (or infinite d) limit. It turns out, however, 
that (|^) has no stationary state in this limit, exhibiting instead a finite-time singularity. By 
adding an extra, physically- motivated term [see (§) below], this deficiency is remedied and 
a stationary state obtained. 

The main properties of the model are: (i) In the limit D — 0+, there is a sharp phase 
transition at a critical coupling Qc from a smooth weak-coupling phase to a rough strong- 
coupling phase; (ii) The strong-coupling phase has a two-peaked height distribution in the 
stationary state, with the height difference depending on g but not on t; (iii) At finite noise 
strength the transition is rounded, but the qualitative features of the weak and strong- 
coupling regimes are very different; (iv) The dynamics of the approach to the stationary 
state is quite different from the dynamics of correlations in the stationary state, contrary to 
the assumptions of conventional scaling theories. 

We shall adopt an approach which parallels the one used in mean-field, infinite range, 
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theories for equilibrium statistical mechanics (e.g. in the Sherrington-Kirkpatrick model of 
spin glasses 0), introduced in this context originally by Luck and Orland ||^. This amounts 
to assuming, once eq. (|1|) is discretized on a lattice, that every site is the nearest neighbor of 
every other site. Therefore the number z of nearest neighbors of each site equals the system 
size. For Euclidean lattices z = 2d, so this approximation is expected to be valid for d ^ 1. 

Let hi{t) = h{xi,t) where Xi is the i^^ lattice point and let y/z be the mesh length of 
the lattice discretization. Let {h)t = lJ2i=ihi(t) be the average height at time t; then 
4>i{t) = hi{t) — {h)t is the local height fluctuation on site i. The discretized form of the 
Laplacian is V'^h\g^ = —(pi, while the nonlinear term becomes {'Vh)'^\^. = ((0^)t + 0^)/2. 
After subtracting the equation dt{h)t = ^{4>'^)t for the mean height, (|I|) becomes 

dt4> = -i^^ + \{(p'-{^')t)/2 + ( (2) 

where the site subscript has been suppressed. Here the noise, in view of our choice of the 
lattice discretization, satisfies (C(i)C(^')) = '^DS{t - 1') with D = Dqz"'^/^. 

That (0) leads to difficulties at large t can be understood by expressing it in the form of 
a Langevin equation, 

d,<j) = -V'i<P)+v, (3) 

(here and below the prime denotes a derivative) and noting that the potential V is unbounded 
from below as — oo. This suggests that there is no stationary state, and that a finite-time 
singularity occurs. The 'self-consistent' time dependence of V which comes from the 
term in (0) make things more complex. Numerical integration of the Langevin equation 
(evolving an ensemble of 2; ^ 1 systems simultaneously to compute confirm the 

conclusion of a finite time singularity. This suggests that the d —>■ 00 limit of the KPZ 
equation is, as it stands, not well defined. 

Since the KPZ equation is usually obtained from a small gradient expansion, it is rea- 
sonable to try to remedy to the failure of eq. (^ by including a further term. In order to 
do this we need to specialize to a particular model. Indeed while the terms appearing in eq. 
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(H) belong to the small gradient expansion of a large number of models, which is the origin 
of the assumed KPZ universality, the first neglected term differs from model to model. For 
example, in the small gradient expansion for directed polymers, the next term in (|I]) would 
be k(V^/i)^. On the contrary, in models for growing interfaces the next term is proportional 
to (VhyV'^h Focusing on the latter models, in the following, we will discuss the equation 

dth = iyV^h + Xlyhf + K{vh)\^h + r]. (4) 

Grouping the first and the third term in the right hand side of eq. @, this reduces to eq. 
with a space dependent coefficient u which is larger on steep portions of the interface (k > 0). 
This is appropriate for restricted solid-on-solid (RSOS) models for growing interfaces, 
where, for entropic reasons, the fluctuations are suppressed faster by the constraint on more 
inclined surface elements [0. The second term also results from the effect the constraint 
has on the growth mechanism. Incoming particles are more likely to stick to the interface, 
causing it to grow, on flatter parts of the surface than on inclined ones . Eq. (|) includes 
both these effects of the constraint to the same (second) order in their gradient expansion 

i- 

Let us now discuss the infinite range limit of eq. for k > 0. Performing the same 
steps as before leads to 

9*0 = fit-il + st)<p + g<p^ - 0' + C- (5) 

Here we have set v = k = 1 via the transformations t ^ t/u and (f) (f)^u/ k. Accordingly 
g = \/{2^/i7k) and D = Dqki'~'^z~^^'^ , where Dq is the strength of the noise in eq. (^. 
Furthermore in this equation we have defined the (time-dependent) parameters 

fit = {<P')t-9{<P')t and St = {<P')t (6) 

which have to be determined self-consistently. The parameter nt enforces the condition 
{4>)t = for all t. The interface velocity is Vt = dt{h)t = gst — Ht- 

Eq. (H) can be easily put in the form (H) with a time-dependent potential 
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vm = '^^t<^ + {l + st)\-g\ + ^ (7) 

which is now bounded from below. 

If we assume the parameters ^ and s approach finite hmits as t — *■ oo, then the stationary 
state is given by 

Pi=oo(0) (xexp[-V;o(0)/^]. (8) 

On the other hand, the values of /Xoo and s^o which enter the potential Voo are obtained from 
eq. (H) from the distribution Poo(0)- A stationary state exists if one can find fioo{g,D) and 
Soo{g,D) which solve self-consistently equations (|^, |^) and (H). 

Let us consider first the limit D ^ 0^. For small g the potential V^o has a single minimum 
at 0^ = 0, giving Soo = = jjoo- For g > ^/S, the potential has a pair of inflection points 
and for g > 2 a second minimum at (p^ = ■^{g + \J g^ — 4) appears. However, for g < 3/\/2 
this minimum has a higher potential than the minimum at = 0, so its occupation is 
exponentially suppressed with respect to that of the other minimum for D ^ 0~^, and Soo 
remains zero (in fact s^o = 0{D) for D small). It is only when, increasing g, V"oo(0^)— Voo(0~) 
becomes of order D that the weight of the + minimum becomes comparable to that of 
the - one. In the limit L> = 0+ this happens at 5'c(0+) = 3/^/2 = 2.12132.... For 
g > gc both minima contribute to (0^)oo so that Soo starts to increase with g even for 
D = 0+. Since both minima must have non-vanishing occupation to ensure (0)oo = 0, s^o 
and /ioo adjust their values such that both minima have the same depth. This implies that 
Voo{(p) = V + ^(0 — 0^)^(0 — (p^y where 0^ are given by V^'oo(0^) = 0. Moreover, in this 
limit Poo(0) = a(5(0 — 0^) + (1 — a)5(0 — 0^). Using (0)oo = and evaluating (0^)oo one 
finds Soo = ^0 0^- The solution to these equations yields (p^ = ^g ± \/|fi'^ ~ | and 

^00(^,0-^) = ^-^, /^oo(^,0+) = -|(^-i) (9) 

for g > gc- The potential at the bottom of both minima is = — These minima are 
separated by a maximum at 0° = (?/3 with Voo{4>°) = -^g^ig^ — 3). 
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The emerging picture for D = is that of a sharp phase transition at Qc- This separates 
physically different weak-coupling {g < g^) and strong-coupling {g > g^) regimes in which 
the static and dynamic properties differ substantially. Note that D = 0+ means we are 
taking the limit D ^ {] after the limit t oo, in order to establish a stationary state. 

Before commenting on the nature of the two phases, it is worth discussing the behavior in 
the critical region at nonzero noise. This is available if one includes the gaussian fluctuations 
in the saddle point calculation of eqs.(|]). If one defines gdD) = 73 (''" ~'~ ^ 
e = g ~ gc, the result can be cast in the scaling form Soo = Df{e/D), for e and D both 
small, where f{x) satisfies 3(1 + / + In a// — 1) = V2x. A similar scaling form holds for 
/ioo = — v^-D — Soo/V^- For D > the transition is smoothed and a crossover between the 
two regimes occurs. We checked by a numerical integration of the Langevin eq. (j^) for ~ 10^ 
realizations (see fig. |lD that the stationary state discussed so far is the one which is actually 
reached by the dynamics. 

Let us discuss our findings in the framework of kinetic roughening. The parameter s has 
the meaning of the square of the interface width s = {{h — (h))^). Its behavior with the 
linear system size L, in the stationary state, defines the roughness exponent x- Soo{L) ~ L^^. 
In our case, we considered from the beginning an infinite system, so a finite value of s^o, both 
in the weak and in the strong coupling phases, implies x = 0. This agrees with the result 
found for directed polymers in random media in the d ^ 00 limit [|l^. However, the nature 
of the stationary state differs dramatically. Indeed eq. (|D predicts, for the strong coupling 
regime, a double peaked distribution in contrast to the highly asymmetric but single peak 
distribution found for directed polymers |Tl|]. The most striking feature of the distribution 



d) is that it predicts the presence of bumps of a definite height 0+ — (f)~ in the surface for 
g > gc- Moreover, the minima of the potential are separated by a barrier AV^ ~ (s'/S)^ 
which is very large for (7 ^ 1. This suggests very large relaxation times r ~ exp{AV/D) 
in the stationary state. The numerical simulation of eq. (^ indeed confirms this conclusion 
(see fig. |I|). The interface velocity, Voo = gs^o — fJ'oo, also differs in the two regimes, being 
0{D) for g < gc and 0{g - gc) for g > gc. 



Let us now discuss the non-stationary dynamics in the strong-couphng regime. Consider 
the Langevin dynamics of an ensemble of a 2; ^ 1 system {(pi} starting from a flat substrate: 
Pt=o{(f)) = Therefore sq = = fiQ. It is convenient to think again in terms of the 

potential (|^). Initially this has two minima, but only the minimum at 00 = is occupied, 
while the one at (p^ ~ g, which is very deep, with Vo(0^) — —g'^/12, is empty. As time 
goes on, the locations of the minima, which are given by V/(0^) = 0, adjust according to 
the running values of St and At the same time, the potential barrier between the two 
minima increases. After an initial transient, the distribution of the ensemble {(pi} can be 
well approximated by Pt{(p) = <ytS{(p — (pt) + (1 ~ C(t)S{(p — (pi). This allows one to transform 
the equations {(p)t = and St = {(p'^)t in the form at = ~(pl /{(pf — (pt) and St = —(pt(pt- 
These relations leave only one independent variable, whose dynamics is given by 

Vt{(P1) - VMT 



da 

—— ~ exp 
dt 



(10) 
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This equation describes the noise activated processes (p^ across the potential barrier. 

It neglects the backward processes (p^ —>■ (p~ . It is therefore valid as long as Vt{(p^) — 
Vt{(p^) ^ D. Since the scale of Vt is of order g'^ ^ this is a good approximation for times 
for which the system is not too close to the stationary state. 

The equation (|10]) cannot be solved exactly. However some progress can be made as- 
suming a <^ 1, and considering the leading behavior in a: (p'^ = g{l — a), 0" ^ —g<y and 
(p° ^ 2ga. To leading order in a -C 1, after some algebra, one finds 



D 2D \tc 



where the characteristic time = (fp') has been defined. Note that in the stationary 
state a^o —>■ 0.1464 . . . as g 00. Therefore the leading order in at can give a reasonable 
approximation even for late times. 

For at ^ tc ^ 1, the saddle point approximation to Pt{4>) becomes questionable. In 
the regime tc ^ a* ^ 1, where we expect our approximations to be reasonable, eq. ([TI1|), to 
leading order in at, predicts a behavior at ~ (Int)^/^. The "interface width" St can easily be 



related to at. Within our approximations {g^ ^ D and at ^ 1) one finds St — D + {(p'^Yat 
to leading order in at- Therefore our previous discussion suggests 

St~(lnt)i/3. (11) 

Figure |I| shows that for = 3 and D = 0.0625 eq. (|TT|) yields an accurate fit of the simulation 
results. 

This regime holds in general for D <^ g'^ (when the saddle point approximation is valid) 
and for intermediate times. When the backward processes 0"*" 0^ becomes important 
(y(0~) — ^(0^) ~ -D), eq. ( p!0D ceases to hold. The system starts to relax to the stationary 
state (see figure p. 

It is important to note that the 'off equilibrium' dynamics has no relation whatsoever with 
the stationary dynamics. In particular we found that the relaxation times to the stationary 
state are orders of magnitude smaller than the ones which characterizes the dynamics in 
the stationary state. In the case illustrated in figure |1|, for example, the fit of the approach 
to the stationary state yields a characteristic time 340, whereas the correlation function 
in the stationary state decays over characteristic times of order exp(AV^/Z}) 8000. The 
scahng relation usually assumed to relate the early stages of growth to the stationary state 
behavior ^ fails dramatically. Finally we observed that the time scales involved in the 
dynamics of the model, even for moderate values of g and D, can well be beyond the reach 
of numerical simulation (already for (7 = 4 and D = 0.125 the time to reach the stationary 
state is of order 10^). 

In conclusion, we have discussed a high- dimensional limit for growth equations. We have 
shown that this limit is plagued by finite time singularities for the KPZ equation as it stands. 
We have modified the model by including a further term in the small gradient expansion of 
the equation for RSOS growth models. Even though this term would be irrelevant, according 
to the usual folklore based on dimensional analysis, we have shown that it allows for a well 
defined and non-trivial high dimensional limit. In particular, in the stationary state the 
resulting surface is characterized by bumps (or depressions, for g < —gc |^) of a definite 
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height in the strong-couphng regime. Although the analysis has been carried out explicitly 
in the limit of vanishing noise, D = 0^, the same qualitative picture obtains for D > 0, with 
the height distribution changing (but now smoothly) from single-peaked to double-peaked 
with increasing g. 

This strong-coupling stationary state differs substantially from the one of directed poly- 



mers in random media. Our results suggest that the mapping discussed in ref. ||T2[ between 
RSOS models and directed polymers may fail in high dimensions. One might wonder whether 
the scenario outlined so far is specific to eq. @ or whether it applies in a more general class 
of models. Note, in this respect, that the sign of k plays a very important role. For k, < 
the infinite dimensional limit is plagued by finite time singularities like those found for the 
KPZ equation. Actually, eq. with k < would result from the small gradient expansion 
of the equation for the Eden model or for ballistic aggregation [0. On the other hand, for 
other generalizations of the KPZ equation (like e.g. the one used in ref. |]13[) in the spirit of 



RSOS models, we found the same scenario as in this work |]T^. These observations raise the 
question of whether a KPZ universality class exists at all in high enough dimensions. It is 
interesting to note that already for = 2 there are discrepancies between the exponents ob- 
tained by numerical simulations for different models |T5|. It is also noteworthy that already 
in d = 2 finite-time singularities in the simulation of the KPZ equation have been observed 
0- 

If the behavior discussed above applies to any large but finite dimension, a second key 
question is what is the dimension above which this behavior sets in. In this respect we 
note that finite height bumps and structures with very large relaxation times have been 
observed in c? = 2 simulations of a regularized KPZ equation. 



We finally point out that the existing numerical results on RSOS models [|T^] in high 
dimensions do not completely rule out the possibility that the behavior we have discussed 
sets in for relatively low dimensions. Indeed in ref. |18| only the early stages of growth were 
investigated and an exponent was extracted under the assumption of dynamic scaling. Our 
model shows that this assumption is far from trivial. Furthermore, in a numerical simulation 



for the interface thickness it might be quite hard to disentangle a power law with a small 
power from a behavior like that of eq. (pJ]). 

We thank T. Blum, M. A. Moore, and T. J. Newman for discussions. This work was 
supported by EPSRC grant GR/H27496. 
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FIGURES 

FIG. 1. Upper curve: Stationary state correlation function C{t) = {(f)(t + to)^(io)) for a system 
with (7 = 3 and D = 1/16. The fit to C{t) = ^exp(— t/r) yields Inr ~ 8.97 to be compared to 
AU/-D = 9. Lower curve: behavior of st, starting from the flat configuration, for the same system. 
The fit was obtained plotting s^2±0. 015 + 0. 017 In f versus hit. 
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